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We analyze the effect of bichromatic microwave irradiation on the magnetoresis- 
tivity of a two dimensional electron system. We follow the model of microwave 
driven Larmor orbits in a regime where two different microwave lights with different 
frequencies are illuminating the sample {w\ and W2) ■ Our calculated results demon- 
strate that now the electronic orbit centers are driven by the superposition of two 
harmonic oscillatory movements with the frequencies of the microwave sources. As a 
result the magnetoresisitivity response presents modulated pulses in the amplitude 
with a frequency of Wl ~ W2 f whereas the main response oscillates with Wl + W ' 2 , 

PACS numbers: 
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In the field of Condensed Matter Physics, very few problems have produced such a 
intense activity, experimental and theoretical, like Microwave Induced Resistivity Oscilla- 
tions (MIRO) 1 ^ and Zero Resistance States (ZRS) 3 ' 4 . From the experimental standpoint, 
remarkable contributions are being published in a continuous basis. Among them we can 
highlight activated temperature dependence in the magnetoresisitivity (p xx ) respons o^ 4 ^ , 
quenching of p xx response at high microwave (MW) intensities 2,6 , absolute negative con- 
ductivity (ANC) and breakdown of ZR S 2 ' 5 ' 6 ' 7 , suppression of MIRO and ZRS by in-plane 
magnetic field^ and the evidence that MIRO and ZRS are notably immune to the sense 
of circular polarization of MW radiation 1 - . Very recently, an experimental achievement 
has joined the group of new contributions. In this case an ultraclean 2DES is subjected 
to bichromatic MW radiation coming from two monochromatic sources with different 
frequencies w\ and wt^. The unexpected result consists in a p xx response which is 
clearly modulated in the oscillations amplitude. This modulation results to be tunable 
by increasing intensity of one of the MW sources keeping the other constant. All these 
experimental evidences establish real challenges for the theoretical models presented to 
dat e 12 ! 18 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 . Considering that all these models are not able to achieve con- 
sensus about the true origin of these striking phenomena, the new experimental results 
can be regarded as crucial tests for theories, for the existing ones, and for the ones to 
come. However some theoretical contributions are already being presented giving expla- 
nation for some of the experimental outcomes. We can stress theoretical proposals for 
temperature and high MW intensity dependence 2 ^ 2 -, absolute negative conductivity^ 3 -^, 
p xx immunity to the sense of circular polarization of MW radiation 25 ' 26 and finally one 
very recent proposal regarding p xx response to bichromatic MW radiation^ 7 . 

In this letter we report a theoretical explanation to the p xx modulated amplitude of 
a 2DES when is illuminated for two different MW radiations with different frequencies 
Wi and u>2- In a recently presented model by the authors 19 , it was demonstrated that 
a 2DES subjected to a perpendicular magnetic field and MW radiation, Larmor orbit 
centers oscillate with the same frequency as the MW field: MW driven Larmor orbits. 
If we now apply two MW fields to the sample the consequence is that electronic orbit 
centers are subjected simultaneously to two oscillatory movements with the frequencies 
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of the MW fields w\ and w 2 . The outcome is the superposition of both oscillations giving 
rise to a oscillatory movement in the orbits center whose amplitude is modulated in the 
way of pulses. Pulses have the well-known frequency of Wl ~ W2 whereas the new oscillatory 
movement goes with Wl + W2 . 

Following the MW driven Larmor orbits model, we first obtain the exact expression of 
the electronic wave vector for a 2DES in a perpendicular magnetic field B and two MW 



sources 



19,28,29. 
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where 7 is a phenomenologically-introduced damping factor for the electronic interaction 
with acoustic phonons^i, e is the electron charge, <p n is the solution for the Schrodinger 
equation of the unforced quantum harmonic oscillator, W\ and w 2 are the MW frequencies, 
Ei and E 2 are the intensities for the MW fields, w c is the cyclotron frequency, X is the 
center of the orbit for the electron motion and x c \ is the classical solution of a forced 
harmonic oscillator driven by two different time dependent forces: 
Xdif) = Xi(t) + x 2 (t) = 



COS W\t + 



eE 2 



Ai cos Wit + A 2 cos w 2 t 



: COS W 2 t 



L is the classical Lagrangian and J m and J n are Bessel functions. If W\ is not very 
different from w 2 and the MW fields intensities are equal, then we can write, 
A\ ~ A 2 = A and therefore: 

x c iit) = A [cos Wit + cos w 2 t] = 2Acos \\{wi — w 2 )t\ cos \_\{wi + w 2 )t\ 
showing that now the oscillatory movement for the Larmor orbits center presents modu- 
lated amplitude with a frequency given by \{wi — w 2 ) whereas the main oscillation goes 
like |(wi+W2)- 

Now we introduce the scattering suffered by the electrons due to charged impurities 
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randomly distributed in the sampl e 1 ? 1 ^ . Following the model described in ref. [19], firstly 
we calculate the electron-charged impurity scattering rate 1/r (being r the scattering 
time). Secondly we find the average effective distance advanced by the electron in every 
scattering jump, that in the case of two MW sources is given by: 

AI = AX° + A\ cos W\T + A2 COSW2T, where AX° is the effective distance advanced 
when there is no MW field present. Again if A\ ~ A2 = A we can write: 
AX MW = AX° + 2Acos - w 2 )r] cos [^(w 1 + w 2 )t] . Finally the longitudinal con- 

ductivity a xx can be calculated: a xx oc J dE AX (fi — //), being fa and ff the corre- 
sponding distribution functions for the initial and final Landau states respectively and 
E energy. To obtain p xx we use the relation p xx = 2 g " 2 ~ ^f 2 -, where a xy ~ ^ and 

^xx'^xy ^xy " *-> 

@xx @xy 

In Fig. 1, we present calculated results for experimental frequencies^ (A7GHz and 
31GHz). In panel a) we represent all the cases (bichromatic and both monochormatic) 
together for comparison. In the rest of panels we represent each calculated response 
separately. According to our model the surprising profile obtained for p xx response is a 
reflex of the amplitude modulated oscillatory movement of the Larmor orbit centers, when 
they are under the influence of both MW fields. In fact in the experimental graphs^ and 
in the calculated Fig 1., at least one modulated pulse can be seen. This pulse can also 
be observed more clearly in Fig. 2, where we represent same situation as in Fig 1., but 
with different frequencies: (57GHz and 38GHz). Our results are in good agreement with 
experiments. 

In Fig. 3, we represent in the bottom graph calculated bichromatic magnetoresistivity 
p xx as a function of B, for experimental frequencies at increasing MW- intensity for the 
case of 47GHz. In the top graph of Fig. 3, we present mathematical functions which 
simulate p xx behavior for monochromatic and bichromatic driving forces, considering one 
of them with increasing magnitude. In both graphs it can be observed that a minimum 
shows up when one of the driving forces is increased. This is because increasing only 
one of the intensities (47GHz in this case), the bichromatic response will tend to be 
similar to the monochromatic one. We have considered that, according to our model 
(see ref. [19]), p xx depends on B like p xx oc B cos(wt) = £> cos (S^) , where S c is a 
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sample dependent scattering term. In the case of two simultaneously MW driving forces 



with similar intensities: p xx oc B cos 



1 (Wx— Wj) 

2 B 



COS 



1 (w 1 +w 2 ) 

2 B 



. Considering the last 



expression we have tested our current theory using the experimental data (Fig. 2 of 

ref. [11]). The aim is to eventually obtain the experimental frequencies used (47GHz 

and 31GHz). In order to do that first we have measured in the experimental figure 

the B~ l periodicity in monochromatic and bichromatic graphs and also the width of 

the corresponding bichromatic pulse. Once obtained this information, and using our 

expressions for the p xx dependence with B we have been able to reach numerical values 

for the pulse and main oscillation frequencies. For the pulse we have obtained a value of 

9.3GHZ, and for the main frequency 38.9GHz. Comparing these values with the ones 

obtained directly from the experimental frequencies, ( 47 ~ 3 ^> GHz = gGHz for the pulse 
and ^+3i)GH Z = 39GHz for the 

main frequency, we can see that the agreement finally 

achieved is quite reasonable. 

In conclusion we have demonstrated that the experimental results 13, regarding the p xx 
modulated response of a 2DES subjected to a B and bichromatic MW can be explained 
in terms of the electronic orbit centers being driven by the superposition of two harmonic 
oscillatory movements with the frequencies of the microwave sources. 

This work was supported by the MCYT (Spain) grant MAT2005-06444, the "Ramon y 
Cajal" program (J.I.). and the EU Human Potential Programme HPRN-CT-2000-00144. 
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Figure 1 caption: Calculated magnetoresistivity p xx as a function of B, for experi- 
mental MW-frequenciesii: AT GHz and 31GHz. a) Calculated MW responses for all the 
frequencies considered, monochromatic and bichromatic MW sources, all together for 
comparison . b) Bicromatic response [47GHz and 31GHz). c) Monochromatic 47GHz. 
d) monochromatic 31GHz. 

Figure 2 caption: Same as Fig. 1 but for MW frequencies 57GHz and 38GHz 

Figure 3 caption: a). Bottom graph: calculated magnetoresistivity p xx as a function of 
B, for experimental frequencies at increasing MW-intensity for the 47GHz MW source. 
Top graph: Mathematical functions to simulate p xx behavior for monochromatic and 
bichromatic driving forces, considering one of them with increasing magnitude. 
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FIG. 1: Calculated magnetoresistivity p xx as a function of B, for experimental MW- 
frequenciesii: 4:7GHz and 2,1GHz. a) Calculated MW responses for all the frequencies consid- 
ered, monochromatic and bichromatic MW sources, all together for comparison . b) Bicromatic 
response (47GHz and 2>\GHz). c) Monochromatic 47GHz. d) monochromatic 31GHz. 
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FIG. 2: Same as Fig. 1 but for MW frequencies 57GHz and 38GHz. 
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FIG. 3: a). Bottom graph: calculated magnetoresistivity p xx as a function of B, for experimental 
frequencies at increasing MW-intensity for the 47GHz MW source. Top graph: Mathematical 
functions to simulate p xx behavior for monochromatic and bichromatic driving forces, consider- 
ing one of them with increasing magnitude. 



